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Garrett, The Jaipur Observatory and its builder, presented by the 
1 author; Greenwich Astrographic Catalogue, vol. ii., 16 charts of 
the Astrographic Chart of the heavens, and Observations of the 
planet Eros, 1900-1901, presented by the Royal Observatory, Green¬ 
wich ; Professor E. S. Holden, Galileo, and other tracts, presented 
by the author; Oxford Astrographic Catalogue, vol. iv,, presented by 
the University Observatory, Oxford; Pennsylvania University 
Publications, Catalogue of 648 double stars discovered by Professor 
Hough, presented by the University; three lantern slides of the 
Corona of 1908 January 3, taken by Professor Campbell, presented 
by the Lick Observatory. 


An Empirical Law of Astronomical Refraction , 

By H. H, Turner, D.Sc., F.R.S., Savilian Professor, 

1. The following investigation was originally undertaken with the 
view of substituting a simple proof of the law of refraction for 
students who could not follow the more elaborate proof involving 
the differential equation. But it was a surprise to find how closely 
the observed refraction could be represented with so rough a sup¬ 
position as that of three, or even two, homogeneous shells of atmo¬ 
sphere ; and the question was suggested whether, in the present state 
of our knowledge, more elaborate hypotheses were really justified. 
If a rough supposition fits the facts, clearly it is no proof of the 
correctness of a more elaborate one that it also fits the facts, 

2. Moreover, suspicion of the correctness of existing hypotheses 
was suggested from another direction. Meteorologists are finding 
that the temperature of the atmosphere does not follow a smooth 
gradient, as is generally assumed in refraction hypotheses; at a 
certain height a wholly unexpected state of things has been found 
to exist. According to M. Teisserenc de Bort there is above 10 or 
12 km,, an “isothermal layer” in which the temperature ceases to 
fall as we ascend; and the conditions are similar over parts of the 
world where the temperatures close to the surface differ widely, It 
is difficult to reconcile these results of observation with the hypo¬ 
theses usually adopted in constructing tables of refraction.* It is, 
however, not intended to examine at present the consequences of 
taking M. Teisserenc de Bort’s work into account;—merely to 
show that it may not be difficult to do so when we have fuller 
information, without, perhaps, dislocating existing refraction tables. 

3. First let us consider an atmosphere of one homogeneous 
spherical layer. Let C be the Earth’s centre; OM a section of its 
surface through C and the star ; LB the section of the boundary of 

* I am indebted to Mr. Saunder for a reference to a paper by Professor 
Bakhuyzen in KoninTc. Akad. van Weten. te Amsterdam, 1907 January 26 ; 
see Nature , 1907 April 4, p. 538; in which the discordance between observation 
and the assumption usually made is demonstrated by a concrete example, 
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the atmosphere ; ABO the direction of the star’s light, refracted 
once only at the boundary B. The star is thus seen in the 
direction OBD, instead of in the direction BA. 

If ABE = <£, -^OBC = \f/, ABOL = f, 

then by the law of refraction 

sin <£> = /& sin ij/ . . . . (i) 

and from the triangle OBC 

sln ^ = bc sm ^ = r+i sin ^ ' ' • (2 ) 

if 2 be the ratio of the height of the atmosphere to the 
Earth’s radius. 




Hence sin cj> = — tL^ sin £ , . . . (3) 

4. The similarity of equations (2) and (3) suggests a geometrical 
construction for the refraction which is obvious enough, but which 
I do not remember to have seen in print. 

With radius i-ti describe the circle MGK round C as centre 

( fi g- 2). 

Then sin OGC 1 jm 

sin GOL — CG ~~ 1 + z 

Hence OGC = <£, just as OBC = ij/, and the refraction is <£ - \f/ = GCB. 

5. If now we take the case of two homogeneous shells of atmo¬ 
sphere, we see that the refraction will be represented by the sum 
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of two angles GjOB^ G 2 CB 2 , with a gap between; the radii CG 1? 
CB 1 , etc. being of constant length, and the angle ZOG being £, 
the apparent zenith distance, and generally, if there are n layers, 
expressing all the angles of incidence and refraction in terms of £, 
we have— 



/*i 


sin ^ = 7 ^ sin ^=^(7 


8in ^ sln ^=^(7^) sin ^ 

I + Z-, . J M-! . 

in ^ = ^ em ^>i=^(7T^) sin ^ 


sm 


JU-n , pi , _ 

sm^= — S m+ 2 = ^j^- )S mt; 


Mi 


• ( 4 > 


sin xj/ n 


Mi 


M«( ^ 4" %n) 


sin £ 


sin <£„ 


Mi 

Mn+it 1 4* 2„) 


sin £ 




The notation is tolerably obvious: the radii of the successive sur¬ 
faces of separation are i+z^ i+z 2 , 1 +z 8 * • • • i+z n -; and the 
refractive indices of the successive homogeneous shells are jm v 
• * *. • /x n . If the n th be the outermost, fx n+z = 1* 
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6. It will be seen that although B x (fig, 3) lies ou the boundary 
of the first shell, B 2 does not lie on the boundary of the second. 


The radius 


CB 2 is 


not 


i+z 2 but — (1 + z 2 ); and if the refractive 

index decreases with the height, this is less than 1 +z 2 . 

Recent meteorological experiments, however, suggest .that there 
may be anomalies in the upper strata. If the refractive index for 
the r th layer should increase with height instead of decreasing, then 
CG r would be greater than CB n and the angle G r CB r would be 
subtractive. 

7. It is of interest to consider whether there are always gaps 
,between the elements. If the angle B 1 CG 2 vanishes, then </* 2 = i/q, 

or /h(i+« 1 )"/ A 8(i+*2) * * • • ( 5 ) 

In this case the refractions of the two assumed layers are 
exactly equivalent, at all zenith distances, to that of a single 

layer with outer radius 

CB 2 
.CG, 


1+2, 


a 

CB 2 = — (i + z 2 ), and of refractive index 
fh 


1 +z 1 ^ 

8. Condition (5) may be regarded as an equation determining 

m 3 , the refractive index of the shell surrounding tie two inner 
Wes j and hence we see that, given any'two homogeneous spherical 
shells of atmosphere with any refractive indices , they may he re¬ 
placed for purposes of astronomical refraction hy a single homo¬ 
geneous shell , provided the refractive index of the medium im¬ 
mediately surrounding them has a certain value , and in that Case 
only . , 

This case must be carefully distinguished from the one in 
which a really homogeneous shell is arbitrarily treated as two by 
an imaginary surface of separation, at which there is no refraction. 
In this latter case one of the angles GjCBj or G 2 CB 2 would be 
zero} not necessarily the gap BjCGg. 

9. Attention is drawn to the-possibility abovo italicised, because 
Jt seemed at first that it might afford the explanation of the curious 

fact that the rough hypothesis of two or three homogeneous shells 
of air can give so gqod an approximation to our observed refractions 
as below. If it is possible to replace two elements by a single one, 
and to continue the process, there is nothing unreasonable in the 
exact representation of a large number of elements by a simple 
equivalent. But the numerical values obtained do not support this 
view, the gaps being wide compared with the thickness of the shells. 

iQ. We now turn to numerical evaluation of the constants. 
The law of refraction is simple enough at moderate zenith distances, 
and we shall be chiefly concerned with its form near the horizon, 
la that case <£, i/r, and £ are all nearly 90°, and it is more con¬ 
venient to deal with their complements, which are small angles. 
Let us therefore write 

£ = 90° — A, 4 > = 90°- v, .\j/ = 9o°-A, 
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/ 

Let a and /3 be the values of A and v for £ = 90° or A = o°. Then 
from equations (2) and (3), 

1 ^ /x . 

COS a = —;—COS p =—:— . . . (6) 

i + z’ ^ i -f z x ' 

and equations (2) and (3) may now be written, 

cos X = cos a cos A, cos v — cos f 3 cos A . . (7) 

In fig. 2, ^OCH = a and ^lOCK = f 3 : ^KCH = a - 0 . The re¬ 
fraction is r = <j> - \j/ = X - v . 

Hi In terms of the method usually adopted, we should find our 
constants from what we could learn of //, and z. From the point of 
view of this note, we are to determine a and (3 empirically from 
observed refractions at two different zenith distances; and the most 
convenient Z.D.’s to select are 45 0 and 90°. We shall confine 
pur attention to mean refractions; and it is convenient to have the 
horizontal refraction exactly 33', for a reason which will become 
clear when we deal with three separate layers. The mean Pulkowa 
refractions, as shown in the second column of Table I., have there¬ 
of been multiplied by a constant. Hence one equation between 
a and (3 is 

a ~ P — 33 = '00960 . . - (8) 

12. At altitude 45 0 the equations may be combined to obtain 
the usual approximate law of tan Z.D. Subtracting one of 
equations (8) from the other, we have 

(cos v — cos X) = (cos f 3 — cos a) cos A , 
or 

• A, — v / n x * A + v 

2 sin-= (cos (3 - cos a) cos A cosec-. 

2 2 

As a first approximation, we can neglect the differences between 

A, v, and A on the right, and put A = 90° - £. If the re- 

2 

fraction be denoted by r, we have for moderatp zenith distances 

o . A - v 

r = 2 sin-. 

2 

. \ r = (cos f 3 - cos a) tan £. 

So that (cos (3 — cos a) is the well-known u coefficient of refraction ” 
at moderate zenith distances: and when £ — 45°, tan£ = 1, and 

cos / 3 - cos a = 55"*7o = •0002700 . * (9) 

Hence 

^ . a — B . a + B 

2 sin- 1 - sm- lL = *0002700 . 

2 2 
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13. Substituting the value of a- ft from equation (9), we find 

i( a + /3) = I ° 3 6 ' 43 i 
and since i( a ~P) — 16 3° 

a =i° 53 ' 13" j8-i° 20'13" . . (10) 

14. We can now estimate the error of the approximation above,. 

where A is substituted for At A—.45** 

log cos v = log cos 4^° - *000118 
log cos A = log cos 45 0 - *000236 
v- 45 ^ 56 ", A - 45 0 = 1 52". 

These differences are found at once from a table of log cosines. 
Hence 

A + v AF ,° v f _ jf 

-= 45 1 2 4 ) 

2 • . 

and the accurate equation to replace equation (9) is 
cos p - cos a = *ooo2 700 x 1*00041 
*(“ + £)= r ° 3 6 ''43 "+i"'9. 

= 1° 3 6 ' 45 " sa y. 

<*=i°53' 15" 20' i S " . . (11) 

15. If a table of logarithms of cosines to seconds of arc is not 
available, we can throw the equation 

cos A = cos a cos A 

into the form 

cos A — cos A 1 — cos a 
cos A + cos X 1 + cos a 

, A — A, A + A , 9 (t 

tan -. tan-= tan J -. 

222 

A+ A 


or 


Let 


tan 


A-A 


= —: then since /v ~ r — = A + 
22 2 


A-A 


PiiPi 4" 2 tan A) = 2 tan 2 -(2 — p x tan A). 


px= - a± Ja 2 + 4b 
4 b 


a+ Ja 2 + 46 


, the upper sign being the 


only one to take; where a = tan A sec 2 -, b — tan 2 - 

2 2 


Pi* 


2(1 - COS a) 


tan A + Jt&n 2 A + sin 2 a 


(12) 
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When we neglect sin a in comparison with tan A, the denominator 
becomes 2 tan A, and 

Pi = ( 1-cos a) cot A. 

Similarly, 

p 2 = (i ~ cos fi) cot A, 

and r — p 1 — p 2 = (cos /? — cos a) cot A as before. 

When sin a is small compared with tan A we have 


and when A =45' 


p, - (1 - cos a) cot A |^i - | • • 0 3) 


= (1 _ COS a) (1 - J sin 2 a) 


a 2 a , 4 


= — — — to the second order in p,. 
26 1 


_/ 3 2 p 4 . 

and p 2 are equal to their tangents. Hence the refraction is 

Pi ~ P2 ~ ( a fi) • —~ ^( a2 • 

Since a — /3 is the horizontal refraction we have 


Similarly, p 2 = 2 —cL • and to" this order of approximation p x 
2 • o 


(h) 


k- 


Refraction at 45 0 
Refraction at Horizon 


= ^[r-^+ /?*)] . (IS) 


16. Let us now consider what happens if we have n concentric 
shells characterised by 

( a l fi\)r~(P“>fi2) * ■* ° 0 '( a n) fin)' 

In each case the horizontal refraction is ( a r - / 3 r ); so that the total 
horizontal refraction is 

( a l ~ fil) + ( a 2 ~ fit) + ( a 3 “ fiz) + • • • + ( a n - fin) = 33- 

How the division of our shells is somewhat arbitrary. We can 
make them all of; equal thickness, for instance; or we can choose 
their relative thicknesses so that they have equal refractions at 
some specified zenith distance. There is a convenience in making 
their horizontal refractions all equal , so that ’■ ‘: 

“l - A = - Pi = a -3 - A)-= • ... = a n~Pn = ^( 33 ') • OO 

Proceeding to the refraction at 45 0 , it is 

2(a,-A)f ttr + A 


(i7) 


and hence substituting for a r - fi n we have 

nk = Z (“""3 /3 ’ r )[ 1 - + & 2 ) ■ ■ 
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where h denotes, as in equation (15), the ratio of refractioh at 45 0 
to that afc the horizon. 

One Shell. 

17^ The numerical results for one shell (with <x=i° 53' 15", 
/?= i°<20 15" as in equation (i;)) are shown in the third column 
of Table I., and the discordances from the Pulkowa refractions in 
the fourth column. It is clear that, though the accordance is close 
as far as Z.D. 6o°, the Hypothesis does not fit the facts at all. The 
discordance at Z.D. 89° is nearly 5', or 20 per cent. 

Table I. 


Calculated Two Shells. 

- A - 


Z.D. 

Pulkowa 

Refraction. 

Calculated 
one Shell. 

0—c. 

Inner. 

Outer. 

-A. 

Sum. 

o-c. 

45 

557 


// 

557 . 

0*0 

9*5 

// 

46*2 , 

. 557 

0*0 

6°. 

. 96*2 


96*2 

0*0 

16*3 

79*9 

96*2 

0*0 

70 

151-9 


152*4 

- 0*5 

26*0 

126*0 

152*0 

: -0*1 : 

r . 3o 

• 305'5 


311*8 

- 6*3 

: 53*5 

253*5 

307*0 

^. I;, 5 

■ 81 

337*5 


346*0 * 

- 8*5 

39*7 

280 

339*7 


• 82 

376*5 


388 

- ii*5 

66 - 

312 

■ 378 , 

2 ' ~ 

83 

424 *8 


442 

- 1 7 

77 

352 ’ 

429 

-4" 

84 

486*1 


512 

- 24 

89 

402 

491 

-5 

8S 

566*2 


606 

-■ 40 

107 

466, 

573 

-7 

86 

674-2 


738 

- 64 

133 

549 

682 

-8 

87 

8257 


935 

-109 

177 ; 

654 ; 

: 831 

-5 

88 

1049*0 

— 

1239 

—-190— 

260- 

'792 “ ' 

1052 

-3 

89 

1397-0 


1678 

- 281 

470 

925 1 

1395 

+ 2 

90 

1980*0 


1980 

0 

990 

990 , 

1980 

0 

In column 

for one shell, 







a : 

= !° 53 ' 

15" 

log cos 

a= 99997643 

1 



(i 

= 1 20 

15 

log cos 

/? = 9-9998817 


In columns for two 

shells, 






■ 

a i 

0 / 

= 0 41 

0" 

log cos 

“i = 9 , 9999 6 9 i 


- 

■-& 

*='?, 24 

3 Q-'- 

•log CQS 

/>’i = 9 -’ 9999 ? 9 °. 

£ 



a 2 

= 2 48 

: 

log cos 

a ?^ 9 ' 999476 o 





= 2 32 

20. 

log cos 

&=- 9 ' 9995735 ,, . 



Two Shells. 

18. With two shells-we have at first sigldr two more disposable 
constants; really only;one. One of them ^is:used up in the 
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assumption that the horizontal refractions of the two shells are 
equal, as in § 16. Our four constants a x , a 2 , f 3 2 are therefore 
to he determined by the equations 

a i ~~ Pl~ a 2 ~ @ 2 ~ 1 • • * • * * • (18) 

2 * = a i±A[ I -|( ai 2 + /3i2 ) 

and one more condition still to be imposed. Inspection of the dis¬ 
cordances for one shell shown in Table I. suggests that we should 
reduce the largest discordance to zero, i.e . should satisfy the 
observed refraction at 89°. This can only be done by “ trial and 
error ” methods. It was found useful to form tables of cos a cos A 
for values of A= i°, 2 0 , 3 0 . . . io°; and for every minute of arc 
of a. By their use, successive approximations gave the values 

o in o nf n 

a x = o 41 o a 2 — 2 40 30 

&=o 24 30 A,=2 32 o 

as a close approximation; the horizontal refraction of 33/ being 
divided into two of i6J'. 

19. We must now examine how far these values satisfy the 
equation (19). They are close enough for substitution in the 
square terms. The following small table facilitates this substi¬ 
tution: the unit is *00001. 



Table II. 
Value of £a 2 . 



1 o'. 

20'. 

3°'- 

40'. 

50'. 

6 o\ 

o° 

1 

I 

2 

3 

6 

10 

1° 

14 

18 

23 

28 

34 

4i 

. 2° 

48 

56 

64 

73 

82 

91 

3° 

101 

112 

124 

136 

149 

162 

4° 

176 

I90 

205 

221 

237 

254 

5° 

i 

271 

289 

3°7 

326 

346 

366 

Equation (19) 

runs 






3° 13' 25 

" = 32' 45" (1- 

•0000 

S) + s° 

4o' 15" 

(1 - ’OOI 



O f 

= 3 13 

0" — 0 

"*i - i4"‘o. 



Hence our right-hand side is defective by 39". Throwing this on 
to the second shell, which is less effective at small altitudes, we 
take as a trial 

Oj = o 41 o a 2 = 2 48 50 

& = o 24 30 fi s = 2 32 20 
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With these values of the constants the refractions in columns 5 and 
6 of Table I. were calculated \ and it will be seen from column 8 
that the error of their sum is small, 

20. The question left outstanding is the choice of slightly 
different constants which shall perhaps satisfy the observations at 
(say) Z.D. 8o°-87° rather better, at the expense of Z.D.’s 88°-9o°, 
where observations are rare and possibly affected by errors of 
several seconds. In Table III, are given the effects of altering a 1 
and by io', and the similar effects of altering a 2 and p 2 by 10'. 
It will be seen that we can, by increasing the constants for one 
layer, and diminishing those for the other by an exactly equal 
amount, leave the moderate Z.D.’s almost unaffected, and at the same 
time introduce sensible modifications into the large Z.D.’s. 

Table III, 


Change for 10' Alteration of Constants . 


Z.D. 

Inner 

Shell. 

Outer 

Shell. 

Difference. 

Error 

from Table I. 

45 

// 

(3) 

// 

( 3 ) 

0*015 

(o"-o) 

6o 

' (5) 

( 5 ) 

0*065 

(0-0) 

70 

( 8 ) 

• (8) 

0*25 

(o-o) 

8o 

(16) 

(i 5 ) 

2*2 

-I'S 

8i . 

( 17 ) 

(16) 

2-9 

- 2 

82 

20 

18 

2 

- 2 

83 

23 

20 

3 

- 2 

84 

26 

22 

4 

-5 

85 

32 

23 

9 

-7 

86 

40 

24 

.16 

-8 

87 

52 

22 

30 

-5 

88 . 

74 

■ 18 

56 

-3 

89 

100 

7 

93 

+ 2 

90 

0 

0 

0 

0 


21. The first five differences in the fourth column of Table III. 
were formed by a special method* as that used for the later ones is 
unsuitable. We may for moderate Z.D.’s simply differentiate the 
equation 

cos X = cos a cos A. 

Thus 

dX sin a cos A 
da (1 — cos 2 a cos 2 A)* 

. ’ _ sin a 

(tan 2 A + sin 2 a)* 

The refraction due to a thin shell between a and a + 8a at altitude A 
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is — . 8a. The change in this for a change Aa in a is . 8a . A a : 
da . da 2 

■where 

d 2 k _ tan 2 A cos a 
da 2 (sin 2 a + tan 2 A)^ ’ 

and when a is small compared with A 

~ cot A [1 - |(i + 3 cot 2 A) sin 2 a]. 

, da 2 


If we have two shells with characteristics cq and a 2 , the difference 
in effect will be proportional to 

J cot A (1 +3 cot 2 A) (sin 2 a 2 - sin^). 

The value of the coefficient £ cot A + § cot 8 A is given in the 
following small table. - 


Table IY. 


Yalne of C = (£cot A + f cot 3 A) and D = o" *007560. 


* Z.D. 

Alt. 

c. 

D. 

0 

45 

O 

45 

2*00 

u 

0*015 

6o 

30 

8*6 

0*065 

70 

20 

32*6 

0*25 

8o 

10 

286 

2*2 

8i 

9 

380 

2- 9 

82 

8 

545 

4 ’ 1 , 

^3 

7 

813 

6*i 

84 

6 

1295 

9-8 

85 

5 

2240 

16*9 

factor for the coefficient C 

in our special 

case is 

10' xarc of 

16J' x (sin 2 2 

i 40 — sm 2 0 

33 ') = o"- 9°756 


The values of D for Z.D. 82°-85° are larger than those found 
directly in Table III., probably because the approximation used in 
Table IY. is ceasing to 'be accurate; but the quantities are of the 
same order of magnitude. 

22. Table III. makes it clear that we cannot sensibly improve 
the accordance at 85° without introducing a much larger error at 
89°. Every 1" of alteration at 85° gives 10" at 89°. Hence we 
may accept the constants used in Table I. as giving'as good a 
representation of the Pulkowa refractions as we can fairly get with 
two shells of ^qual horizontal refraction. r; • , . , 
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Three Shells. 

23. We divide the horizontal refraction equally into three of x i'. 
It is an obvious convenience to have whole minutes of arc, and for 
this reason the horizontal refraction was taken as exactly 33' in the 
first instance. We might have taken 36', which divides by 2, 3, 
and 4 ; but some work had already been done with 33' before it 
was found that two shells would give such a good approximation. 

After some few trials, the following values were found to be 
near the truth :— 

O / // O / .// O t It 

a x ^=o 29 o a 2 =i 24 o a g = 3 14 p 

A = ° 18 0 & =I 1 3 o /? s = 3 3 o 

i( a i + A) = ° 33 3 ° K a 2 + ^2) = I 18 3 ° i ( a 3 + A)) = 3 8 3 ° 

Hence equation (17) becomes in this case 

4 0 50' 8" = 23' 30" (1 - *00003)+ i B 18' 30" (1 - *00034)+ 3 0 8' 30" 
(1 - *00191) 

' = 4 0 5 o # 30 #/ -o / '-o 4 -i w - 5 -ai / '-6 

showing that the adjustment has been properly made. 


Table V. 

Three Homogeneous Shells. 


Z.D. 

Inmost 

Shell. 

Middle 

Shell. 

Outmost 

Shell. 

Sum = C. 

o-c. 

Difference for 10 '. 

Inmost. Middle. Outmost. 

45 

4‘47 

15*09 

36'14 

5570 

// 

0*00 

n 

195 

u 

1*85 

i-8 5 

6o 

7-67 

26*16 

62 *42 

96-25 

0*00 

3'56 

3*18 

3*i8 

70 

12*21 

41*34 

98-27 

151-82 

+ 0*1 

5 ' 3 ° 

5*05 

5*05 

8o 

25’3 

85*0 

i960 

3 ° 6'3 

-0*8 

11-3 

10*2 

9*1 

8i 

28*3 

94*4 

215*8 

33 8 ‘5 

- 1 *o 

12*2 

II*I 

10*2 

82 

317 

106*0 

2397 

377 - 4 

-0*9 

I 3'9 

12*5 

10*9 

83 

36*2 

120*7 

269 *o 

425-9 

-1 *i 

16*2 

14*4 

ii*6 

84 

42*4 

140*5 

305*5 

488*4 

-2*3 

18 *5 

16*1 

ii*8 

85 

51*0 

167*0 

350*5 

568*5 

- 2*3 

22*3 

18*6 

12*5 

86 

63*5 

205*5 

407 *o 

676*0 

-i-8 

277 

21*8 

13*2 

87 

00 

o' 

264*0 

477*5 

826*5 

-o*8 

39*0 

28*5 

n*5 

88 

i 25'5 

362*0 

556-5 

1044*0 

+ 5*o 

• 52*3 

28*6 

8*2 

89 

236-5 

525*o 

629*0 

i39o*5 

+ 6*5 

83*2 

19*5 

2*3 

90 

660 - o 

660 *0 

660*0 

1980*0 

0*0 

o*o 

0*0 

0*0 


24. The results of this hypothesis are shown in Table V., and 
from the seventh column, 0—C, it will be seen that there is 
fair accordance with the Pulkowa refractions. This accordance 
might be improved; by use of the differences given in the last 
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three columns we could form linear equations to find small correc¬ 
tions to the mean radii of the three shells. But the improvement 
would be of the nature of a compromise, and would involve— 

(a) A study of the accuracy of our knowledge of refractions at 
small altitudes. We should have to settle, for instance, what 
error at 89° Z.D. can be made to obtain an improvement of 1" at 
Z.D. 8o°, 

(b) A more accurate computation of the differences for io\ 
Those given above were found in the way made clear by the 
following example;— 

1 st 2nd 
„ diff. diff. 

log cos 18' log COS IO° = 9*9933455 = log cos (io° + 16*0) u 

a 2 5'3 „ 

log cos 29'log COS IO P = 9*9933361 = log cos (io° + 41 *3) I2*4 

o " 377 

log cos 40' log cos io° = 9*9933221 = log cos (io° + 79*0) 

The first two lines are formed to find the refraction due to the 
inmost layer (with constants 18' and 29') for altitude io°, the 
result being 25"*3 as given in the 2nd column. By adding the 
3rd line we find that the refraction of a shell with constants .29' 
and 40' (each 11' greater than the former) is 37^-7 ; so that the 
difference for 11' is i2 ,, *4; and for 10' would be ii"* 3, which is 
given in the 7th column. But in forming differences of differences 
slight inaccuracies are multiplied, and inspection of the figures 
near the top of the last three columns shows that (since the numbers 
are nearly equal for all three shells) we are really concerned with 
differences of these again. 

For the present, no attempt will be made to improve the 
approximation, and we turn to some other points. 

25. It is of some interest to deduce the heights and refractive 
indices of the shells above found empirically, 

One Shell. 

log (1 +z) = log sec a = 0*0002357 = log 1*000544 

log /x = log cos /3 — log cos a = o*oooii74 = log 1*000270. 

The “ height of the homogeneous atmosphere ” is thus indicated as 
4000 x *000544 miles, or 2*18 miles, which is a good deal smaller 
than that assigned by total pressure. 

This is not surprising, for in replacing a series of successive 
bounding surfaces by a single surface, we should expect the equiva¬ 
lent single surface to lie in the midst of the constituents. We can 
imagine, for instance, a uniform graduation which would bring the 
equivalent surface to the mean height, i.e. would give us a £C height 
of the homogeneous atmosphere ” just half that assigned by calcu¬ 
lating the total pressure; and it would not be surprising to find 
that the equivalent surface should be lower than this, 
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Two Shells, 

Equations (4) of § 5 may be put in the form 

Log. 

log (1 + zf = log + log sec oq = o *0000309 
1°8 Ihhh = l°g ( l + z i) + log cos Pi — 0*0000199 
log (1 + z 2 ) = log /ij/x 2 + log sec a 2 = 0-0005439 
log /V/Xj = log (1 + 2 2 ) + log cos ( 3 2 = 0-0001174 


509 


Number, 

1*000071 

1*000046 

1*001253 

1*000270 


If we put //, 3 = 1, then jxj = 1*000270, jjl 2 = 1*000224. The heights 
of the bounding surfaces in miles are ^ = 0*28, ^ = 5*00. 


Three Shells. 


The method of formation of the quantities has been made clear 
in the last case. 

Log, Number, 

log (1 +z 1 ) = o + log sec <^ = 0*0000154 1*000035 

For log fjbJ/ji 2 , add log cos f 3 v 0*0000094 1*000022 

„ log(i+2 2 ) „ log sec a 2 , 0*0001391 1*000321 

>. l°g p-ihh » log COS 0-0000412 1-000097 

„ log(i+z 8 ) „ log sec ag, 0-0007331 1-001689 

„ log „ log COS /? 3 , 0*0001175 1*000270 

If we put /x 4 = 1, we have /x 1 = 1*000270, as in equation (9). 
This gives a check on the work. 

Further, fx 2=1*000248, ^3=1*000173. The heights of the 
bounding surfaces are, in miles, 

^ = 0*14, z 2 = i *28, 23 = 6*76. 


It does not seem likely that these figures can have any physical 
interpretation. They merely emphasise the fact that a serious part 
of the refraction takes place near the Earth, and not high in the 
atmosphere. By comparing the columns for the separate shells 
shown in Tables I. and V. it will be seen how rapid is the increase 
in importance of the inner shells as we approach the horizon. In 
Table V., for instance, the two inner shells together only give half 
the effect of the outmost at moderate Z.D.’s; and they do not pro¬ 
duce so great an effect as the outmost until we reach Z.D. 88°, when 
they surpass and ultimately double it in the remaining 2 0 or 3 0 . 
This suggests that any influence of meteorological phenomena on 
astronomical refractions must be sought at large Z.D.’s, which is in 
accordance with experience. But discussion of such points in detail 
cannot be undertaken at present. 


Summary, 

§§ 1-2. Introductory. 

§§ 3-6. A graphical representation of the refraction of a homo¬ 
geneous spherical shell of atmosphere and of several such shells. 
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§§ 7-8. Condition that two ot more such shells should be replace¬ 
able by an equivalent single shell. 

§§ 10-15. Empirical determination of constants for a single 
shell from observed refractions at Z.D. 45 an( l ^h® horizon. 

§ 16. Extension of formulae to two or more shells. 

§ 17* Numerical results for one shell. The errors at 8o°, 85°, 
and 89° are respectively - 6", - 40", - 281". The hypothesis can¬ 
not be said to fit the facts at all beyond Z.D. 60 or 70 at most. 
See Table I. 

§§ 18-22. Numerical results for two shells. See Table I. The 
errors at 8o°, 85°, and 89° are reduced to - i"*5, - 7", and + 2''—a 
very fair accordance for so obviously rough an assumption. 

§§ 23-24 Numerical results for three shells : see Table V. The 
improvement is not great, much less marked than before; and it 
seems clear that the principal step was taken in passing from one 
shell to two. 

§ 25. Some numerical values for heights, etc. 


The Perturbations of Halley's Comet in the Past . Fourth Paper. 
The period 760 to 1066. By P. H. Cowell, M.A., F.R.S., and 
A. C. D. Crommelin, B.A. 

We have again to acknowledge the kind assistance of Dr. Smart 
and Messrs. F. R. Cripps and Thos. Wright in these calculations. 

In the last paper, M.N., lxviii. 5, p. 378, we found March 27 
for the approximate date of perihelion passage in 1066, and 
44"*686 for the value of n at that date. For the preceding passage 
we used the date given by Hind (989 September 12) for the 
purpose of computing the perturbations, and the following results 
indicate that this date is correct within a few days :— 


Revolution 989-1066. 


Planet. 

Limits of u. 

JtZn. 

Jdcr. 

J*r. 

Venus 

O- 30° 

+ 'Oi9 

// 

+ 532” 

> > 

330-360 

- '0035 

... 

... 

Earth 

0 

i 

00 

0 

+ *007 

... 

+ 196 

>> 

330-360 

- ‘013 

... 

... 

Jupiter 

0- 90 

+ '4535 

“ 53 i 

+12807 

* 7 

90-270 

- *2780 

-203 

“ 2905 

9 9 

270-360 

-1 -2940 

+ 708 

+ 22 

Saturn 

0- 90 

- -0560 

- 2 3 

- 1578 

99 

90-270 

+ *0754 

- 6 

+ 1979 

7 9 

270-360 

- '2334 

+ 7 

+ 47 

Uranus 

0-360 

+‘ *0030 

... 

+ 55 


Sums 

- 1 *3200 

- 48 

+ 11155 
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